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ABSTRACT 

An extended XMM-Newton observation of the Seyfert 1 galaxy NGC 4051 has revealed 
a rich absorption line spectrum indicating the presence of a photoionised outflow with 
a wide range of velocities and ionisation parameter. At low continuum fluxes an emis- 
sion line spectrum is well defined with both narrow and broad emission components 
of several abundant metal ions. The absorption line velocity structure and a broad 
correlation of velocity with ionisation parameter are consistent with an outflow sce- 
nario where a highly ionised, high velocity wind, perhaps launched during intermittent 
super-Eddington accretion, runs into the interstellar medium or previous ejecta, losing 
much of its kinetic energy in the resultant strong shock. We explore the possibility 
that a quasi-constant soft X-ray component may be evidence of this post-shock cool- 
ing. This revised view of AGN outflows is consistent with multiple minor Eddington 
accretion episodes creating a momentum-driven feedback linking black hole and host 
galaxy growth. 

Key words: galaxies: active - galaxies: Seyfert: general - galaxies: individual: NGC 
4051 - X-ray: galaxies 



1 INTRODUCTION 

NGC 4051 is a bright, narrow line Seyfert 1 galaxy in the 
Ursa Major cluster, lying at a TuUy-Fisher distance of 15.2 
Mpc (Russell 2003), with a heliocentric velocity of 753 km 
s~^ (Verheijen 2001). The rapid and large amplitude vari- 
ability (Lawrence et al. 1985,1987) was strong early evidence 
that the powerful X-radiation found to be a common prop- 
erty of Seyfert galaxies by Ariel 5 (Cooke et al. 1976, Elvis et 
al. 1978) and Uhuru (Tananbaum et al. 1978) observations 
resulted from accretion onto a central supermassive black 
hole. Subsequent monitoring with RXTE revealed a differ- 
ent aspect of NGC 4051, where it occasionally lapses into 
extended periods of low and quasi-constant X-ray emission 
(Lamer et al. 2003). More recently, high resolution X-ray 
spectra provided by Chandra and XMM-Newton have shown 
that NGC 4051 also exhibits a strong ionised outflow. 

An early Chandra HETG observation resolved two X- 
ray absorption line systems, with outflowing velocities of 
photoionised gas of ~600 and ~2300 km s~^, while contem- 
poraneous HST spectra of CIV, NV and SilV found sev- 
eral absorption systems with velocities from ~30 km s~^ to 
^650 km s~^ (CoUinge et al.2001). Of particular interest in 
the context of the present analysis, where we find a clear 
correlation of velocity and ionisation parameter, the higher 
velocity X-ray component had no counterpart in the UV 



spectra. The Chandra data also showed an unresolved Fe K 
emission line at ~6.41 keV (FWHM <2800 km s"^). 

Observations of NGC 4051 with XMM-Newton in 2001 
and 2002 coincided with periods of relatively high and low 
X-ray flux, offering an opportunity to further explore the 
complexity of its X-ray spectrum. Pounds et al. (2004), here- 
after Po04, found the hard X-ray band to be dominated 
by reflection from cold matter, which could also explain a 
non-varying, narrow Fe K fluorescent line. A soft X-ray nar- 
row emission line spectrum evident at low continuum fluxes, 
with observed wavelengths consistent with the NGC 4051 
rest frame, indicated an extended ionised emission region, 
while a dominant absorption line spectrum in the high flux 
RGS observation of 2001 showed an ionised outflow with a 
line of sight velocity of ~500 km s~^. An Fe K absorption 
line in the EPIC spectrum indicated the presence of a more 
highly ionised outflow component with a velocity some 10 
or 30 times higher, depending on the ionisation state being 
FeXXVI or FeXXV. 

Most recently, Steenbrugge et al. (2009) have reported 
the analysis of Chandra LETG observations in 2001 and 
2003, finding evidence for a more complex ionised outflow, 
with low ionisation components (log^~0.1 and log^^O.5-0.9) 
outflowing at v~200-330 km s~^ , a more highly ionised com- 
ponent of log^'^2 at v~600 km s^^ and a high ionisation 
component (logf~3) at v~4600 km s~^. Those authors find 
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no evidence for recombination in 3 of the 4 ionisation com- 
ponents up to 20 ks after a sudden, factor-of-5 drop in X-ray 
flux, indicating their location at a radial distance ^ 7 x 10^® 
cm for the higher velocity gas and ^ 9 x 10^^ cm for the 
lowest velocity absorption in OV. We note that these mini- 
mum radii conflict with those from a similar time variability 
analysis of Krongold et al. (2007), and have the incidental 
bonus of removing the conceptual problem in Krongold et 
al. where the local escape velocity exceeded the measured 
outflow value. 

Evidence for much higher velocity X-ray outflows 
(Chartas et al. 2002, Pounds et al. 2003, 2006; Reeves et 
al. 2003, Cappi 2006, Tombesi et al. 2010) has been con- 
fined to the very highly ionised matter (log^~3.5-4) most 
readily detected in the Fe K band. The high velocities and 
high column densities in these Fe K band observations are 
thought to offer the best prospect that such energetic flows 
represent a significant feedback mccharnsm to constrain the 
continued growth of a black hole and star formation in the 
host galaxy, although to date only in the case of the bright 
QSO PG1211+143 has the requisite wide angle flow and 
large covering factor been directly observed (Pounds and 
Reeves 2009). In the case of NGC 4051 in contrast, Kron- 
gold et al. (2007) concluded that the low velocity g£is has a 
small covering factor and hence involves a relatively insignif- 
icant mass and energy rate. While Stccnbruggc ct al. (2009) 
do not confirm the high density and small radial distance of 
the iorused outflow, they also derive a low mass and energy 
rate from the soft X-ray spectrum of NGC 4051. 

In this paper we report an analysis of a new XMM- 
Newton study of NGC 4051 with substantially greater sen- 
sitivity than hitherto, finding a rich absorption line spec- 
trum revealing a velocity-structured outflow covering a wide 
range in velocity and ionisation parameter. We re-examine 
the suggestion in Po04 that both high and low velocity ab- 
sorption lines arise at different stages of a mass-conserved 
outflow. Our new assessment considers the new data in the 
context of a high velocity, highly ionised wind being slowed 
by interaction with the interstellar medium, losing much of 
its mechanical energy in the resulting shock, while retain- 
ing momentum to push low ionisation gas out into the host 
galaxy (King 2010). We note that the same ideas would ap- 
ply to an outflow colliding with earlier, slower moving ejecta. 

The wider importance of such shocked outflows, per- 
haps launched during an Eddington accretion episode (King 
and Pounds 2003), lies in the possibility that the accumu- 
lated thrust from multiple episodes - rather than the outflow 
energy - would eventually drive gas from the bulge, thereby 
limiting further star formation and black hole growth. Such 
a momentum-driven feedback mechanism has been shown 
by King (2003,2005) to reproduce the observed correlation 
of black hole and galaxy masses (e.g. Ferrarese and Merritt 
2000, Gebhardt et al. 2000, Haring and Rix 2004)). 



2 OBSERVATIONS AND DATA ANALYSIS 

NGC 4051 was observed by XMM-Newton on 15 orbits be- 
tween 2009 May 3 and .June 15. Here we use mainly the soft 
X-ray data from the Reflection Grating Spectrometer/RGS 
(den Herder et al. 2001), while also checking for evidence of 
a more highly ionised outflow in the EPIC pn (Striider et 



al.2001) and MOS (Turner et al. 2001) spectra. Excluding 
high background data near the end of each orbit the net 
exposures available for spectral fitting were typically 40 ks 
per orbit, yielding an overall exposure of ~600 ks (pn) and 
~1.2 Ms (combined MOS and RGS), factors of ~6 and ~12 
higher than in the XMM-Newton observations of 2001 and 
2002. In the present paper we use spectra integrated over 
whole orbits, and note the mean orbit flux levels cover a 
similar range to the 'bright' and 'faint' observations of 2001 
and 2002. Full details on the timing, flux levels and X-ray 
light curves for each revolution are included in an accompa- 
nying paper on the 'Rapid X-ray variability of NGC 4051' 
by Vaughan et al. (2010). 

Figure 1 contrasts the RGS spectrum at representative 
high (revs 1727-1730) and low (revs 1725,36,39) flux levels 
during the 2009 observation. The distinction found previ- 
ously (Po04) is again evident, with a dominant absorption 
line spectrum seen against the stronger continuum when 
bright, and an emission line spectrum emerging more clearly 
when the continuum is faint. The most obvious absorption 
lines are associated with H- and He-like ions of C, N, O and 
Ne, and multiple velocity components are seen for several 
resonance transitions. The emission spectrum is also better 
defined than in any earlier observation of NGC 4051, bene- 
fiting from the long exposures of the 2009 observation. 

Reference to figure 1 also shows a broad trough in the 
RGS spectrum between ~15.5-16.8 A, present in both high 
and low flux data, which we attribute to a UTA from Fe-L 
absorption in low ionisation matter. Inclusion of this UTA 
casts doubt on the reality of a relativistic emission line of 
OVIII Lyman-a, previously reported for NGC 4051 by Ogle 
et al. (2004) and Steenbrugge et al. (2009). 

In analysing the 2001 XMM-Newton absorption spec- 
trum of NGC 4051, Po04 first subtracted line emission de- 
termined in the low flux spectrum of 2002. However, in the 
2009 XMM-Newton data the low flux spectrum still shows 
residual absorption. We therefore adopted a different ap- 
proach, by first modelling the emission spectrum in the low 
flux data and retaining the measured emission parameters 
in fitting the high flux absorption spectrum. 

A representative 'highflux' spectrum was obtained by 
integrating RGS 1 and RGS 2 data from the 4 orbital revo- 
lutions 1727-1730, while broad resonance emission lines and 
radiative recombination continua (RRC) were determined 
from a representative 'lowflux' spectrum (revs 1725, 1736 
and 1739). The 'highflux' continuum was first modelled by 
a power law (r^3.5), with broad emission features added 
as observed in the low flux spectrum. Individual absorption 
lines were then fitted with negative Gaussians, with wave- 
length, line width and amplitude as free parameters. 

We discuss the emission line spectrum of NGC 4051 in 
more detail in a separate paper (Pounds and Vaughan 2011; 
hereafter Paper II). For the present purposes, in order to 
quantify the absorption line spectrum in the highflux state, 
we simply represent the broad line emission component of 
each principal resonance line with a positive Gaussian, with 
parameters determined in the corresponding lowflux spec- 
trum. 
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Figure 1. Fluxed RGS spectra of NGC 4051 summed over 4 relatively bright and 3 faint continuum levels illustrating the change from 
a complex absorption line spectrum to a dominant emission line spectrum as the continuum flux falls 



3 A COMPLEX ABSORPTION LINE 
SPECTRUM 

3.1 OVIII Lyman-a 

Figure 2 (top panel) illustrates the lowflux spectrum in the 
region of the OVIII Lyman-a line (rest wavelength 18.97A). 
A broad emission line is clearly resolved, with substantial 
absorption near the line core. Fitting a positive Gaussian 
to the broad emission component finds a rest wavelength of 
18.96± 0.02 A, Icr line width of 90± 9 mA and amplitude of 
3.7 X 10~* photons cm^^ s^^ A^^. The median wavelength 
indicates a velocity blueshift of ~750±200 km s~^ while, 
allowing for the RGS 1 resolution {a— 28 mA at 19A), the 
line width corresponds to --3100±400 km s"^ FWHM, if 
interpreted purely as velocity broadening. 

Figure 2 (lower panel) shows the same spectral region 
near the OVIII Lyman-Q line, with a complex of absorption 
lines superimposed on the continuum plus broad line emis- 
sion. A sweep across the 18-20 A waveband with a negative 
Gaussian of width initially set to the RGS resolution, found 
3 strong and 3 weaker absorptions lines, located at wave- 
lengths labelled (a) - (f) in the figure. However, it is ap- 
parent that the absorption structure is only approximately 
matched by Gaussian lines, with the broad opacity between 
-18.5-18.8 A suggesting a spread of velocities. 

Identifying the absorption components (a) - (d) with 
OVIII Lyman-a yields line-of-sight velocities of -780±70 km 
s"S -4100±85 km s-\ -5600±60 km s"^ and -9000±300 km 
s~^. We note that the probable blend of component (c) with 
the OVII ls-3p line (rest wavelength 18.63 A) renders the 
parameters of that component uncertain. However, velocity 
components a, b and c are all independently supported by 
their detection in OVIII Lyman- /3 (Table 2). 

The weaker absorption lines labelled (e) and (f) in figure 

2 are most likely due inner shell transitions in OVI (rest 
wavelengths 19.18 A and 19.34 A; Holczer et al. 2009), and 
we note a further possible satellite line in OV (19.97 A), all 

3 low ionisation lines yielding a low outflow velocity. 

In summary, we find evidence for 3, or possibly 4 blue- 
shifted velocities in OVIII Lyman-Q, with a clear separation 
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Figure 2. (top) A broad emission component of OVIII Lyman-a 
(rest wavelength shown by the vertical line) in the lowflux spec- 
trum of NGC 4051 together with narrow emission and absorp- 
tion components, (lower) Gaussian fitting to absorption lines ob- 
served in the highflux continuum spectrum plus broad emission 
line. Gaussians labelled a, b, c and d are attributed to outflow 
components of increasing velocity in OVIII Lyman-a. Gaussians 
e and f are identified with low velocity absorption in OVI 
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between the component at ^^780 km s ^ and components at 
-4100 and ~5600 km s"\ 

3.2 OVII ls-2p triplet 

The lowflux spectrum in the region of the OVII ls-2p triplet 
also shows a strong and broad emission feature (figure 3, top 
panel), together with a narrow emission line readily identi- 
fied with the forbidden transition (rest wavelength 22.101 
A). Fitting a Gaussian to the narrow line yields a wave- 
length (adjusted for the redshift of NGC 4051) of 22.092± 
0.003 A, line width (t=28± 3 mA and amplitude 9.5 x 10"'' 
photons cm~^ s~^ A~^. We find the OVII forbidden line to 
be marginally blue-shifted, with a velocity of -125±40 km 
s~^, and to be unresolved by RGSl, the 1 a upper limit 
corresponding to a velocity width of ~250 km s~^ FWHM. 

A single Gaussian fit to the broad emission from the 
OVII ls-2p triplet finds a rest wavelength of 21.8±0.1 A, 
line width (t=305±35 mA and amplitude 2.2 x 10~* pho- 
tons cm~^ s~^ A~^. The intrinsic width, corresponding to a 
velocity width of ~9500 km s~^ (FWHM), is clearly much 
larger than that found for OVIII Lyman-Q. However, a more 
detailed analysis of the emission spectrum in Paper II finds 
broad emission components for both the OVII resonance 
(2I.6O2A) and intercombination (21.807 A) lines, with com- 
ponent line widths comparable to that found here for OVIII 
Lyman-Q. For the present purposes of quantifying the high- 
fiux absorption spectrum we retain the single emission pro- 
file, shown added to the best-fit continuum in figure 3 (lower 
panel) . 

Six absorption lines are detected in the region of the 
OVII triplet, all but that near 22.75 A being of width consis- 
tent with the RGS resolution. Those labelled (a), (b), (c) are 
identified as blue-shifted components of the OVII ls-2p res- 
onance line, the individual blue shifts yielding line-of-sight 
outflow velocities of -440±60 km s"\ -4080±140 km s"\ 
and -5810±120 km s~^, respectively. While these values are 
consistent with the 3 stronger components seen in OVIII 
Lyman-Q, the high velocity absorption is much weaker in 
OVII. We note the strong low velocity component in OVII 
is probably a blend of core absorption in the broad emission 
line (confirmed as self absorption in Paper II) and absorp- 
tion in the highflux continuum. 

Deep absorption lines are also identified with inner shell 
transitions in OVI (d), OV (e) and OIV (f), confirming the 
presence of substantial low ionisation gas in the outflow from 
NGC 4051. The deduced velocities for these low ionisation 
lines are comparable to the lowest velocity component in 
OVII. Detailed values are listed in Table 2, apart from line 
(f), where a probable blend of OIV lines at 22.74 A and 
22.77 A (Holczer et al. 2009) makes a velocity measurement 
uncertain. 

3.3 NVI ls-2p and CVI Lyman-/3 

Given the complexity of the emission and absorption in the 
region of the OVII triplet, it is interesting to examine the 
RGS data for other He-like ions. That for CV unfortunately 
lies beyond the range of the RGS, while the spectrum in the 
vicinity of the NelX triplet is complicated by the presence 
of strong Fe-L lines. That makes NVI the best candidate for 
direct comparison with OVII. 



o - 

CO 

E 
o 




20 21 22 _ 20 

observed wavelength (A) 

— 1 — 1 — 1 — 1 — I — 1 — 1 — 1 — 1 — I — 1 — 1 — 1 — 1 — f 




c b a d e f 



20 21 22 23 

observed wavelength (A) 

Figure 3. (top) A broad emission line in the lowflux spectrum 
at the OVII triplet is a blend of components from the resonance 
and intercombination lines. An unresolved emission component 
is identified with the forbidden line in the OVII triplet, (lower) 
Gaussian fitting to absorption lines in the highflux spectrum in 
the wavelength region of the OVII triplet. Those labelled a, b 
and c are identifled as blue-shifted components of the OVII ls-2p 
resonance transition (rest wavelength 21.60 A), while components 
d, e and f are identifled with the low velocity flow in OVI, OV, 
and OIV, respectively 



Figure 4 shows the highflux data in the region of the 
NVI triplet. As for OVII, a broad emission line has been 
added to the power law continuum as a baseline. The for- 
bidden line (rest wavelength 29.534 A) is again unresolved 
and marginally blue-shifted. 

The strongest absorption lines, labelled (a) and (b), are 
readily identified with the low velocity outfiow component 
seen in the ls-2p resonance line of NVI (rest wavelength 
28.787 A) and CVI Lyman-/? (rest wavelength 28.466 A), 
respectively. The deduced velocities of -325±50 km s~^ and 
-260±40 km s"'^, the lowest found in the analysis so far, 
extend the observed link of velocity with ionisation energy. 

We note that any higher velocity component of NVI, 
corresponding to that observed near 4000 km s~^ in OVIII 
and OVII, would be blended with the absorption in CVI. 
Possible interpretations of the absorption lines (c), (d), (e) 
include lower ionisation stages of N, but the identifications 
and rest wavelengths are not sufficiently secure to derive 
useful fiow velocities. 
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Figure 4. Absorption lines in the liighfiux spectrum reveal low 
outflow velocity components of the NVI ls-2p resonance line (a) 
and CVI Lyman-/3 (b). The vertical lines again correspond to the 
rest wavelengths of these transitions 



3.4 Lyman-a of CVI, NVII and NeX 

Outflow velocities found in OVIII Lyman-a at ~5600 and 
~9000 km are considerably higher than any previously 
reported from NGC 4051. However, the statistical signifi- 
cance of the highest component in OVIII is marginal, un- 
derlining the need to look for evidence of a similar veloc- 
ity structure in other H-like resonance transitions available 
within the RGS waveband. Those are of CVI, NVII and 
NeX. In modelling the highflux absorption spectrum, a shal- 
low emission line of width a=90 mA is added to the best-fit 
highflux continuum in each case. 

Figure 5 (top panel) shows the RGS spectrum in the 
waveband 32 - 36 A, which covers the region of CVI Lyman- 
a, and the ls-3p transition of CV. The same 4 velocity com- 
ponents detected in OVIII Lyman-a may be seen as ab- 
sorption lines (a), (b), (c) and (d) in CVI Lyman-a (rest 
wavelength 33.74 A), with values here of -490±35 km s~^, - 
4000±70 km s"\ -5900±90 km s'S and -8920±180 km s"\ 
Identifying the absorption line (e) with a low outflow ve- 
locity in CV ls-3p (rest wavelength 34.97 A), of -250±60 
km s~^, further extends the trend where the lowest outflow 
velocities have the greatest opacity in the lowest ionisation 
stages. 

Figure 5 (mid panel) shows absorption lines in the spec- 
tral band encompassing the NVII Lyman-a line (rest wave- 
length 24.781 A), and the ls-3p resonance line of NVI (24.90 
A). A Gaussian sweep finds 5 absorption lines. Those la- 
belled (a), (b), (c) and (d), yield outflow velocities of - 
495±50 km s'^, -4000±80 km s-\ -5920±135 km s-\ and 
-8600±120 km s~^, consistent with values seen in OVIII and 
CVI. For NVII Lyman-a, the deduced velocities are . Com- 
ponent (e) corresponds to an outflow at -325±50 km s~^ if 
correctly attributed to NVI. The positive Gaussian at ~25.3 
A is an approximate match to the strong RRC of CVI. 

The flnal Lyman-a line accessible in the RGS wave- 
band is that of NeX, having the highest ionisation potential 
and shortest wavelength (rest wavelength 12.13 A). Figure 
5 (lower panel) shows the highflux absorption spectrum. Al- 
though the RGS suffers from decreasing sensitivity in this 
short wavelength region, several absorption lines are seen. 
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Figure 5. (top) Absorption line spectrum in the wavelength re- 
gion of CVI Lyman-a and the ls-3p transition of CV. Lines la- 
belled a, b, c and d are identified with velocity components in 
CVI Lyman-« consistent with velocities found in OVIII. Line e is 
identified with a low velocity outflow in CV. (middle) Absorption 
lines in the wavelength region of NVII Lyman-a and the ls-3p 
transition of NVI. Lines labelled a, b, c and d are identified with 
velocity components in NVII Lyman-a similar to those found in 
OVIII and CVI. Line e is identified with a low velocity compo- 
nents of NVI ls-3p. (lower) Absorption lines in the wavelength 
region of NeX Lyman-a. Lines labelled a and b correspond to the 
lowest and highest velocity outflow components found in OVIII, 
CVI and NVII. Line c corresponds to an intermediate velocity in 
NelX, while lines d, e and f are probably identifled with FeXX 
or FeXIX, but the uncertain identification limits their value as a 
velocity indicator 



6 K. A. Pounds et al. 



Line (a) lies just shortward of NeX Lyman-a (rest wave- 
length 12.13 A), and is most probably due to a low velocity 
component of -620±150 km s~^. If also attributed to NeX, 
line (b) corresponds to a much higher outflow component at 
-8230±130 km s^^, while intermediate velocity components 
at ~4000 and ~6000 km s"^ are not detected in NeX. We 
suggest the most likely identification of line (c) is with NelX 
ls-3p (rest wavelength 11.55 A), at a velocity of -5000±350 
km s~^ It is interesting that for NeX only the highest and 
lowest velocity components of those seen in OVIII, NVII 
and CVI are detected, a further indication of the trend of 
increasing outflow velocity with ionisation energy of the par- 
ent ion. 

Absorption lines (d), (e) and (f) are probably identified 
with FeXX or FeXIX, but the uncertain individual transi- 
tions makes their present value as a velocity indicator not 
very useful. 

In summary, the velocity structure is similar in all 4 
Lyman-Q lines visible in the RGS spectrum, all showing ve- 
locity components up to ^8000-9000 km s~^, in addition 
to a well-separated low velocity component. The highest ve- 
locity of ~8500 km s"'^ is notably strongest in the highest 
energy Ne ion, extending the trend of velocity being corre- 
lated with ionisation parameter at least in the high velocity 
gas. We note, in contrast, that the low velocity absorption 
is seen in all resonance lines across the RGS spectrum. 

In Section 6 we suggest that this difference, also re- 
flected in XSTAR modelling, is due to a separate origin of 
a component of the low velocity absorption arising from self 
absorption in a limb-brightened shell, rather than in the di- 
rect AGN continuum. 



3.5 A velocity- and ionisation-structured outflow 

Table 2 brings together all the outflow velocities assigned to 
a significant absorption component where the parent ion can 
be reliably identified. Also listed in Table 2 is the ionisation 
parameter at which each ion stage would have a maximum 
concentration in a photoionised gas illuminated by the AGN 
continuum flux (Kallman et al. 1996). Several conclusions 
can be drawn from the data in Table 2, which are also shown 
graphically in figure 8. 

First, several discrete outflow velocities are indicated, 
each across a range of difi^erent ions, with a clear distinction 
between a 'low velocity' group (~250-800 km s^^) and a 
higher velocity group, with components at ~4000, ~6000 
and ~8500 km s~^. Second, there is a strong correlation of 
velocity with ionisation parameter, absorption in the lower 
velocity components being strongest in the lower ionisation 
lines, and the higher velocity components being noticeably 
stronger for the higher ionisation lines. Moreover, only the 
lowest velocity flow is evident in the lowest ionisation gas. 

To ascertain whether the trend of increasing outflow 
velocity with ionisation parameter continues to still higher 
values, we look for evidence of absorption lines in the EPIC 
data, covering the K-shell energies of the heavier abundant 
elements Mg, Si, S. Ar and Fe. 



4 ABSORPTION LINES IN THE EPIC DATA 

4.1 Intermediate mass ions of Mg to Ar 

Examination of the EPIC spectra reveals several absorption 
lines, detected in both pn and MOS data, lying close to K- 
band resonance transitions in Mg, Si, S, Ar and Fe. 

Figure 6 shows a section of the highflux (revsl727-30) 
pn and MOS spectra in the energy range from 1-4 keV. Al- 
though the plotted ratio of data to a simple power law con- 
tinuum fit are rather noisy (with deviations near 1.8 keV 
and 2.3 KeV that are probably due to imperfect modelling 
of the Si K and Au M edges in the instrument response 
function), a consistent set of absorption lines is seen for the 
4 abundant elements with K-shell resonance transitions in 
this energy range. 

Identifying the features, observed in both pn and MOS 
data near 1.50, 2.05, 2.68 and 3.42 keV, with Lyman-a ab- 
sorption lines of Mg, Si, S and Ar, respectively, and taking 
a weighted mean of Gaussian fits to the MOS and pn data, 
yields outflow velocities of -7100±1500 km s"\ -7100±2000 
km s-\ -7700±1600 km s'^ and -9700±1800 km s'^ km 
s^^. Within the substantially larger uncertainties these val- 
ues are consistent with the higher velocity components in 
the RGS data, extending their detection over a still wider 
range of ionisation parameter. 

The EPIC velocities are added to those from the RGS 
analysis in Table 2. 



4.2 The Fe K region 

Narrow absorption lines are also seen in the Fe K region 
of the EPIC spectra. Correct modelling of the continuum 
is especially important in the neighbourhood of the Fe K 
lines as spectral curvature below ~6.4 keV will have an as- 
sociated Fe-K absorption edge if attributed to reflection or 
partial covering. Cold reflection was shown in Po04 to dom- 
inate the hard lowflux spectrum of NGC 4051 in 2002, a 
conclusion consistent with the 2009 observation. It seems 
a reasonable assumption that a similar underlying compo- 
nent will contribute to the spectral curvature - and the Fe 
K emission line - in the present highflux data. 

Figure 7 shows a ratio plot of the EPIC pn highflux 
data to a continuum fltted at 2-10 keV with a power law, plus 
cold reflection flxed at the value (R'^l) determined from the 
lowflux data. Gaussian fitting finds the Fe K emission line 
can be modelled with narrow and broad components, with 4 
possible absorption lines detected to higher energy. Details 
of Gaussian fits to the Fe K spectral structure are listed in 
Table 1. Of particular interest here are the absorption lines, 
where alternative conclusions can be drawn regarding out- 
flow velocities, depending on the line identifications. Though 
not unambiguous, an appealing interpretation is where 2 ve- 
locity groups are found, with the lower value aligning with 
the higher velocities seen in the RGS data. A still higher 
velocity, close to the value v~0.1c which may be a signature 
of high velocity winds (Tombesi et al. 2010), is also implied, 
though this depends on the significance of the 2 weaker ab- 
sorption features observed at ~7.5 and ~7.7 keV. We return 
to the likelihood of such a high velocity wind for NGC 4051 
in Section 6.2. 

Summarising the EPIC absorption line data, we again 
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Figure 6. Negative Gaussian fits to (top) the highflux pn data 
and (lower) tlie corresponding MOS data, plotted in both cases 
against a power law continuum. Four absorption features common 
to both plots are consistent with an outflow velocity in the range 
6000-10000 km s~^ when identified with the Lyman-a line of - 
from left to right - Mg, Si, S and Ar. Deviations, seen particularly 
in the MOS data near 1.8 keV and 2.3 keV, are attributed to 
imperfect modelling of the Si K- and gold M-absorption edges in 
the detector response and mirror reflectivity 



Table 1. Summary of Gaussian line fits (adjusted to the NGC 
4051 rest frame) to the Fe K spectral structure shown in figure 7 



line 


energy (keV) 


width (eV) 


ident 


velocity (km s 


1 


6.38±0.01 


60 (f) 


Fe-K 




2 


6.4±0.1 


380±63 


Fe-K 




3 


6.81±0.02 


60 (f) 


FeXXV 


-4900±750 


4 


7.10±0.01 


60 (f) 


FeXXVI 


-5800±450 


5 


7.49±0.03 


60 (f) 


FeXXV 


-31500±2000 


6 


7.72±0.03 


60 (f) 


FeXXVI 


-29500±2000 



emphasise that the deduced outflow velocities depend on the 
correct identification of less well resolved absorption features 
(than for the RGS spectra). Identifying the absorption fea- 
tures in Mg, Si, S and Ar with the relevant Lyman-a line 
is the most conservative, in the sense of yielding the lowest 
velocities, showing the intermediate velocity outflow seen 
in the RGS spectra extends to a still higher ionisation pa- 




observed energy (keV) 



Figure 7. Gaussian fitting to spectral structure in the Fe K band 
of the highfiux pn data, plotted against a power law plus refiection 
fit to the continuum, models the Fe K emission line with broad 
narrow and broad components, while 4 possible absorption lines 
are detected to higher energy. Details are listed in Table 1 

rameter, log^~3. Interpreting the Fe K absorption is more 
ambiguous and we retain both the intermediate and high 
velocity values in Table 2. 



5 MODELLING WITH XSTAR 

To further quantify the photoionised absorption in the com- 
plex outflow in NGC 4051 we have modelled the highflux 
RGS spectrum with XSTAR. We again fitted the contin- 
uum with a power law, adding a negative Gaussian line to 
represent the broad UTA evident in figure 1 near 16 A. As 
in the Gaussian fitting, broad emission features (line plus 
RRC) were added to the continuum. 

A sequence of photoionised absorbers represented by 
grid 18 from the XSTAR library (Kallman et al. 1996) was 
then added to model the absorption line spectrum. Grid 18 
includes a turbulent velocity of 100 km s^^, with an ionising 
power law spectrum similar to that observed for NGC 4051 
irradiating a gas taken initially to be of solar abundance. The 
ionisation parameter, column density and velocity (output 
as a modified redshift) are the primary free parameters of 
each absorber. 

We began by modelling the RGS 1 data over the wave- 
band 17-24 A, dominated by absorption lines of OIV, V, 
VI, VII and VIII and allowing a fit largely independent of 
relative abundances. A visually acceptable fit was obtained 
with 2 photoionised components expressing the low and in- 
termediate velocities. The fit was then extended to 36 A, 
to include the resonance lines of NVII,NVI and CVI, and 
with both RGSl and RGS2 data, and finally over the full 
waveband 10-36 A, covering the higher energy K-shell res- 
onance transitions of Ne, and a potential complex of Fe-L 
lines. For the full band fits the abundances were allowed to 
vary, although tied for the same element across the separate 
ionised components. 

Five photoionised absorbers were required to obtain a 
satisfactory spectral fit over the whole waveband, with a sta- 
tistical improvement from x^/d.o.f. — 6788/4204 to x^/d.o.f. 
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Table 2. Absorption lines identified in the 2009 highflux XMM-Newton spectrum of NGC 4051. The rest wavelength (energy) of each 
absorption line is in Angstroms (keV) and measured line-of-sight outflow velocity components vl, v2, v3, v4, v5 are in km s~^. Log^ (erg 
cm s~^) is the ionisation parameter at which each ion stage would have a maximum concentration in a photoionised gas illuminated by 
the AGN continuum flux. EPIC values are a mean of the pn and MOS flts for Mg, Si, S and Ar while only pn values are shown for Fe. 
Only statistical errors are included 
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u.o 


NVI ls-3p 


24.90 


325±50 


_ 


_ 


_ 


_ 


0.45 


NVII Lya 


24.78 


495±50 


4000±80 


5920±135 


8600±120 




1.1 


OV ls-2p 


22.37 


310±60 










-0.25 


OVI ls-2p 


22.01 


440±75 










0.2 


OVII ls-2p 


21.60 


440±60 


4080±140 


5810±120 






0.9 


OVI ls-2p 


19.34 


850±200 










0.2 


OVI ls-2p 


19.18 


660±200 










0.2 


OVIII Lya 


18.97 


780±70 


4100±85 


5600±60 


9000±300 




1.55 


OVIII Ly/3 


16.01 


500±130 


3700±130 


6300±160 






1.55 


NelX ls-3p 


11.55 






5000±350 






1.4 


NeX Lya 


12.13 


620±150 






8230±130 




2.0 


MgXII Lya 


(1.47) 








7100±1500 




2.25 


SiXlV Lya 


(2.00) 








7100±2000 




2.65 


SXVl Lya 


(2.62) 








7700±1600 




2.9 


ArXVllI Lya 


(3.32) 








9700±1800 




3.2 


FeXXV ls-2p 


(6.70) 






4900±750 




31500±2000 


3.3 


FeXXVI Lya 


(6.96) 






5800±450 




29500±2000 


3.9 



= 5543/4180. The abundances of the most important ele- 
ments (relative to oxygen=l) were C: 0.53±0.16, N:2.4±0.5, 
Ne: 2.6±2.0 and Fe:1.30±0.26. 

The parameters of this multi-absorber fit are hsted in 
Table 3 in the order in which they were added. Components 
1 and 2 evidently represent the high velocity continuum ab- 
sorption observed most strongly in OVIII and higher level 
ions. 

Component 3 represents the low velocity and lower ion 
stage contirmum absorption, as well as the low velocity self- 
absorption in the broad line emission (figure 10). Component 
4 reflects the existence of low velocity absorption across a 
wide range of ionisation parameter. We interpret that in a 
more detailed study of the broad line emission (Paper II) 
where the emission and strong self-absorption arise from 
a limb-brightened shell, the velocities being limited by the 
high- inclination to the line of sight (figure 10). 

Component 5 is intriguing in indicating a small rcdshift. 
We suggest that the red-shifted absorption, which impor- 
tantly can also be seen in the data, arises from the deceller- 
ating outflow on the far side of the AGN, the low ionisation 
parameter being consistent with an origin in the higher den- 
sity gas accumulating ahead of the contact discontinuity (see 
Discussion) . 

Components 1-4 are added in figure 8 to a plot of out- 
flow velocities obtained from the Gaussian line fitting, where 
each ion is located at the ionisation parameter where it 
would have maximum abundance. We have not included the 
uncertainties on each XSTAR component as the fit should 
be considered as representative of a continuous velocity and 
ionisation distribution. Nevertheless, as discussed in Section 
6, we believe the near linear progression from components 



Table 3. Parameters of the photoionised outflow fitted to the 
RGS data. Column density is in H atoms cm~^ and ionisation 
parameter in erg cm 



Comp 


logj 




velocity (km s ^ 




1 


2.97±0.06 


1.4±0.2xl022 


-5880±60 


484 


2 


2.52±0.07 


1.85±0.4xl02i 


-3850±60 


178 


3 


1.43±0.22 


l.l±0.5xl020 


-530±90 


394 


4 


2.77±0.15 


2.9±1.4xl02i 


-400±50 


64 


5 


0.32±0.24 


1.0±0.2xl020 


-|-120±60 


123 



1 to 3 is an important demonstration of a mass conserved 
outflow. 

A 5th asterisk in figure 8 represents the putative pre- 
shock high velocity flow, with velocity v~0.1c and ionisation 
parameter, adjusted to the current sub-Eddington luminos- 
ity, of log^~3.7. 



6 DISCUSSION 

A striking feature of the absorption spectrum in NGC 4051 
is the very wide range of velocities and ion stages observed. 
A second feature is the correlation of velocity and ionisation 
parameter in the continuum absorption. 

Figure 9 visualises that correlation, with contrasting ve- 
locity profiles of the opacity in OVIII and OVII, for the 
high flux revs 1722 and 1724. The upper plot, centred at 
the rest wavelength of OVIII Lyman-a, shows the onset of 
absorption at ~7500 km s~^, with the opacity increasing to 
~6000-5000 km s~^, whereafter it decreases again to disap- 
pear at ~3500 km , re-emerging strongly below ~1500 
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km s~^. The corresponding velocity profile in the OVII Is- 
2p resonance line shows only weak high velocity opacity at 
^4500-3000 km s~^. In contrast the low velocity absorption 
is notably stronger in the lower energy ion 

Figure 8 brings together the velocity data from Gaus- 
sian fitting to the RGS and EPIC absorption spectra, plot- 
ting each well-defined velocity against the optimum ionisa- 
tion parameter for the parent ion. The strongest absorption 
components at each velocity generally lie to the right hand 

- higher ionisation - side of the plot, consistent with a broad 
correlation of velocity and ionisation parameter, while the 
highest velocities are only detected in the higher energy par- 
ent ions, and the low velocity group are seen most strongly 

- or only - in low ionisation matter. 

The correlation of velocity and ionisation parameter is 
further illustrated by the results from XSTAR modelling of 
the highflux absorption spectrum, represented by asterisks 
in figure 8, with components 1, 2 and 3 following a clear 
linear trend in velocity and ionisation parameter . Figure 8 
also indicates that part of the low velocity absorption does 
not follow the general trend with ionisation, with XSTAR 
component 4 supporting the evidence from Gaussian fitting 
that low velocity absorption exists at all ionisation stages 
covered by the RGS data. Figure 10 shows schematically 
how separate absorption spectra are associated with the di- 
rect continuum and with self absorption in the post shcck 
gas. More detailed consideration of a low velocity absorp- 
tion structure, visible across all ion stages in the RGS data 
(and represented by components 3, 4 and 5 in the XSTAR 
modelling), is deferred to Paper II where it is interpreted as 
self-absorption in a limb-brightened shell. 

Meanwhile, an overview of figure 8 suggests 3 broad ve- 
locity regimes in the outflow in line-of-sight to NGC 4051, 
including the higher velocity component in the Fe K absorp- 
tion. On this overview, the individual velocity components 
picked out by Gaussian fitting at ~4000 km s~^, ^6000 km 
s~^ and ^^8500 km s~^, but appearing more like a broad 
trough in the strongest absorption, in OVIII Lyman-a (fig- 
ure 9), may represent density variations within a continu- 
ous fiow. Further evidence for such density variations, which 
might represent a residual shell-structure linked to the inter- 
mittent nature of the initial fast wind, can be seen in velocity 
profiles taken at different flux levels for different ions. 

In the remainder of this paper we refer to the low 
(<1000 km s"^), intermediate (~3000-9000 km s~^) and 
high velocity (^30000 km s~^) regions indicated in flgure 
8 and explore the possibility that they represent different 
stages in a shocked outflow. The implications of such an 
interpretation are doubly signiflcant. The concept of a slow- 
ing and cooling/recombining ionised outflow is in contrast to 
most current ideas for the radial acceleration of AGN winds, 
while efficient post-shock cooling would mean that the me- 
chanical energy in a fast outflow may not be the primary 
mechanism for AGN feedback. 

6.1 Outflows and feedback in AGN 

Comparison with previous observations of NGC 4051 by 
XMM-Newton and Chandra suggest that the low velocity ab- 
sorption is persistent, while that at higher velocities is prob- 
ably variable, at least on a timescale of years. The caveat 
on variability is that most of the earlier observations were 
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Figure 8. Outflow velocities derived from the Gaussian fitting 
plotted against the optimum ionisation parameter for each par- 
ent Ion stage. Also shown by asterisks are the parameters of the 
4 photoionised absorbers derived from XSTAR modelling of the 
RGS absorption spectra, together with a velocity/ high ionisation 
point to represent the putative pre-shock wind 



of lower sensitivity, while the relatively long 2001 XMM- 
Newton observation does - on a re-examination - show evi- 
dence of blue shifted absorption corresponding to an outfiow 
in the range ~4000-6000 km s"\ The new XMM-Newton 
spectra are most closely approached by the Chandra LETG 
spectra reported by Steenbrugge (2009), who find outfiow 
velocities of -^200, ^^600 and ~4600 km s~^, modelled by 4 
ionised absorbers with ionisation parameters ranging from 

ioge~o-3. 

In discussing the 2001 RGS spectrum of NGC 4051 Po04 
noted that, with the simple assumption of conservation of 
mass in a radial outfiow, an extended region of slow mov- 
ing, low ionisation gas might be a continuation of the high 
velocity, high ionisation flow seen in absorption in the Fe 
K band. On this picture much of the mechanical energy in 
the initial outfiow would have been lost before reaching the 
lower ionisation stage, and Po04 speculated that this might 
be due to internal shocks occurring in the high velocity gas. 

We now take up that idea again but in the context of a 
shock interaction with slower moving matter, either the local 
ISM or previous ejecta. To justify that approach we note 
that the broad correlation of velocity and ionisation state in 
the absorption spectra for NGC 4051 is a clear signature of 
a cooling, decellerating and recombining outflow. 

The mass rate in a radial outfiow of velocity v, and par- 
ticle density n at radius r is Afout = ^irhnr^ .v.nip, where b is 
the fractional coUimation angle and rUp is the proton mass. 
Mass conservation in the fiow requires the product n.r^.v to 
be constant. As n.r^ = Lion/^, if the ionising radiation is 
unchanged (or changes very little), for example over a dis- 
tance small compared with r, then conservation of mass in 
a radial outflow yields a linear correlation of velocity and 
ionisation parameter. We would expect that to be the case 
in a post-shock cooling shell. The important implication of 
a decellerating radial outflow is that the mechanical energy 
in the flow would be substantially reduced as the fiow is 
slowed. 

King (2010) has recently examined a relevant scenario. 
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Figure 9. Contrasting velocity profiles of the opacity in OVIII 
and OVII, for the high flux revs 1722 and 1724, illustrate the 
strong correlation of velocity and ionisation parameter consistent 
with a decellerating post-shock flow. The upper plot, centred at 
the rest wavelength of OVIII Lyman-a, shows the onset of absorp- 
tion at ~7500 km s~^, with opacity increasing to ^GOOO-SOOO km 
s~^, whereafter it decreases again to disappear at ~3500 km s~^, 
re-emerging strongly below ~1500 km s~^. The corresponding 
velocity profile in the OVII ls-2p resonance line shows significant 
high velocity opacity only at '^4500—3000 km s~^. In contrast the 
low velocity absorption is notably stronger in the lower energy ion 



where a highly ionised, high velocity wind drives into the 
interstellar medium, losing much of its energy by efScient 
cooling of the shocked gas. Such a scenario would have ma- 
jor implications for studies of AGN feedback based on X-ray 
absorption spectra. Until now. X-ray observations have been 
used in attempts to show that fast, ionised outflows can pro- 
vide the link between the growth of a SMBH and its host 
galaxy, by the integrated mechanical energy in the fast flow 
(eg Pounds and Reeves 2009). However, as recently pointed 
out (eg King 2010), if such an energetic wind persists while 
the black hole doubles its mass (the Salpeter time), the cou- 
pling of wind energy to galactic baryons may have to be 
inefficient to allow massive bulges to grow to the values ob- 
served, effective feedback instead being enabled by the total 
momentum of the flow, an alternative that has been shown 
to yield the observed M - a relationship for nearby active 
galaxies (King 2003, 2005). 

In the following section the new XMM-Newton observa- 



tions of NGC 4051 are examined in the context of a shocked 
wind, where the intermediate velocity /intermediate ionisa- 
tion outflow corresponds to the immediate post-shock gas 
and the low velocity/low ionisation absorption to matter 
building up ahead of the contact discontinuity. 

6.2 Comparing the NGC 4051 data with a 
shocked wind model 

In the King (2010) shocked wind model a high velocity 
ionised outflow collides with the ISM of the host galaxy, 
resulting in a strong shock. The gas density increases by a 
factor ~ 4 at the shock front, and the velocity drops by 
the same factor. Beyond this (reverse, adiabatic) shock, the 
flow is further compressed in a relatively thin, cooling re- 
gion, while the velocity slows to low values. Strong Compton 
cooling by the AGN radiation implies a fairly rapid transi- 
tion between the immediate post-shock regime and the much 
slower and compressed state near the contact discontinu- 
ity. Beyond the contact discontinuity a further low velocity, 
low ionisation component will build up as the interstellar 
medium is swept up by an outer (forward) shock. 

Comparison of that scenario with the new NGC 4051 
data assumes an intermittent highly ionised wind, with typ- 
ical values of v~0.1c and log^~4 (Tombesi et al. 2010), has 
collided with the ISM or previous slower moving ejecta, with 
the density increase rendering the immediate post-shock gas, 
v~0.025c (~7500 km s'^), visible as resonance absorption 
at log^<3. 

The shocked wind cools and slows after passing through 
the inner shock, yielding the broad absorption trough seen 
in the velocity plot for OVIII Lyman-a line (figure 9). We 
speculate that the observed gap between the intermediate 
and low velocity absorbing material arises from the falling 
column density of OVIII, re-emerging at low velocities as the 
column density builds up ahead of the contact discontinuity. 

It is interesting to consider how the observations of 
NGC 4051 may be used to estimate the parameters of an 
earlier, perhaps intermittent, Eddington episode whose ef- 
fects we may now be observing. While direct evidence for a 
high velocity wind is marginal in the 2009 data, as it was in 
the 2001 XMM-Newton observation, Tombesi et al. report 
the detection of an outflow at ~0.13c in the 2002 obser- 
vation. Taken together, those reports indicate the present 
fast outflow in NGC 4051 is intermittent, and we note that 
with the unusually high column density required to detect 
a blue-shifted Fe K absorption line in a low redshift source 
such as NGC 4051, radial expansion would rapidly render a 
transient outflow undetectable with EPIC. 

In what follows, we explore the effects of such an in- 
termittent fast outflow as it interacts with the ISM or 
slower moving ejecta, as an explanation of the velocity- and 
ionisation-structured outflow phases observed in the present 
observation of NGC 4051. 



6.3 The intermediate ionisation/intermediate 
velocity gas 

We identify the immediate post-shock outflow in NGC 4051 
with the onset of substantial opacity at a velocity v~7 x 10* 
cm s~^ (flgure 9, top panel). XSTAR modelling (Table 3) 
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finds an ionisation parameter log^^S and column density of 
Nif~1.4 X 10^^ cm~^ to represent this flow phase. Both ve- 
locity and ionisation parameter are consistent with the fac- 
tor ~4 change from the putative 0.1c wind expected across 
a strong shock. 

The incident ionising luminosity of ~8 x 10''^ erg s~^, 
for an average-flux EPIC spectrum (rev 1729), together with 
the fitted ionisation parameter and measured velocity, give 
nr'^v ~6 X 10'*^ s^^ for the post-shock flow. Assuming an 
angular coUimation b=0.3, the post-shock flow mass rate is 
then Mout ~4x 10^* gm s"^ (~0.06 Mq yr"^). Interestingly, 
that rate is close to the Eddington accretion rate for NGC 
4051 assuming an accretion efficiency of 0.1. 

The corresponding momentum rate of the intermediate 
velocity gas is then ~3 x 10^^ (cgs), with mechanical energy 
(0.5Mout.«^) of ~10*^ erg s~^ That mechanical energy rate 
is ~0.3% of the Eddington luminosity for a black hole mass 
of ~ 1.7 x 10*^ Mq (Denney et al. 2009), a factor ~30 less than 
the value w/cljEdd predicted by a simple continuum driving 
model (King and Pounds 2003) , roughly consistent with the 
expected velocity-linked loss at the strong shock. 

As the post-shock gcis cools we can identify second and 
third XSTAR components on the velocity/ionisation plot 
(figure 8), with v~3850 km s~^ and log^~2.55, and v~550 
km s~^ and log^~1.56, following reasonably closely the ex- 
pectation of a linear correlation of velocity and ionisation 
parameter in a post-shock flow. 

The detection of several strong RRC and broad res- 
onance line emission is indicative of a strongly rccombin- 
ing stage in the flow, and the observed OVII RRC flux of 
~4 X 10~® photons cm~^s^^ provides a measure of the in- 
termediate velocity flow. A temperature of ~5 eV, from 
the width of the RRC, implies a recombination rate for 
OVIII of -10"" cm^ s-i (Verner and Ferland 1996). With 
that value and assuming 30 percent of recombinations from 
the majority OVIII ion direct to the ground state, we de- 
duce an emission measure from the OVII RRC flux of order 
~2 x lO'^^cm"^, for a TuUy-Fisher distance of 15.2 Mpc. 
We note, in passing, that this could be a minimum measure 
as any much higher temperature RRC component would be 
difficult to resolve. 

Importantly, the particle density in the intermediate 
flow region is constrained by evidence of a narrowing of the 
broad emission component in OVIII in rev 1739, indicating 
a change in the ionisation state following a 3-4 day interval 
of unusually low continuum flux level. This is discussed in 
more detail in Paperll. In the following estimates of the flow 
properties we assume a recombination timescale for OVIII of 
~4 days, corresponding to a particle density n~5 x 10®cm~* 
for the intermediate velocity flow. 

The OVII RRC emission measure then corresponds to 
an emitting volume of ~8 x lO^^cm^ and - taking the column 
density of component 2 in the XSTAR modelling as a mea^ 
sure of the intermediate velocity absorber - gives a radial 
thickness of order Ar~4 x 10^"' cm. 

At a mid-phase velocity of —4000 km s~^, the shocked 
gas would traverse this cooling region in a time t— lO'^ s. 
That would suggest the intermediate velocity absorption 
could exist in the absence of a fast outflow for a simi- 
lar timescale, while the considerably shorter recombination 
timescale could explain structure in the intermediate veloc- 




Figure 10. Sketch showing the origin of separate absorption 
spectra, in the continuum by line-of-sight to the AGN and by 
self-absorption in the soft X-ray emission from a limb brightened 
shell 

ity flow, in turn perhaps reflecting the intermittent nature 
of the fast wind. 

Prom the above estimates of particle density, shell thick- 
ness and emission volume (acknowledging that these are only 
a crude measure across a likely strong radial gradient), we 
find a shell radius r— 7 x 10^^ cm, for b=0.3. 

We now compare the above parameters of the interme- 
diate velocity/ionisation post-shock gas with values for the 
low ionisation/low velocity gas accumulating ahead of the 
contact discontinuity, to provide an order-of-magnitude es- 
timate of the duration and history of the putative Eddington 
accretion episode in NGC 4051. 

6.4 Constraining the low ionisation/low velocity 
gas 

We take component 3 from the XSTAR photoionisation 
modelling to represent the low ionisation/low velocity gas 
accumulating ahead of the contact discontinuity (CD) and 
seen in continuum absorption in the highfiux spectrum. That 
component has an ionisation parameter of log^~1.43, ab- 
sorbing column density of ~10^°cm~^, and outflow velocity 
of —530 km s~^. With a mean ionising luminosity ( > 1 keV) 
of 8 x 10*^ erg s~^ , we obtain n.r^— 3 x 10*° for the low veloc- 
ity/low ionisation flow gas. With the particle density scaling 
from the velocity difference - and reflected in the lower ioni- 
sation parameter, we assume n~5 x 10*' cm"^, giving a radial 
depth of the pre-CD shefl Jr - 2 x 10^^ cm"^. 

The higher assumed density of the pre-CD shell would 
appear to conflict with the persistence of the low veloc- 
ity absorption, particular in a higher level line such as 
OVIII Lyman-aZp/ia. However, as we report in Paper II, self- 
absorption of the broad emission line can fully account for 
the low velocity opacity when the continuum is weak. 

As the radius of this low ionisation shell must exceed the 
estimated thickness of the intermediate flow region, we note 
that light travel time delays would further limit variability 
in the low velocity opacity, consistent with the lower limit 
of r>9 X 10^^ cm obtained by Steenbrugge et al. (2009) 
from a lack of variability in the low velocity component in 
the Chandra absorption spectrum. We assume r— 10^* cm 
below. 

With that geometry the mass of accumulated low ioni- 
sation gas, seen in absorption, is ~7 x 10^^ gm (~0.3 M©). 
Comparison with the mass rate in the post shock flow in- 
dicates an accumulation time of ~6 years. For a velocity of 
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530 km s~^, the mechanical energy in this component of the 
low velocity gas is then ~10**ergs, indicating that ~99.5% 
of the mechanical energy has been lost in the post-shock 
cooling. 

While the total flow momentum will be conserved 
through the shock, we expect a major fraction of the ini- 
tial ram pressure is converted to gas pressure at the contact 
discontinuity. Integrating the immediate post-shock outflow 
momentum over 6 years totals ^10*^ (cgs). In comparison 
the accumulated ram pressure in the post-shock, low ionisa- 
tion gas is ~3.7x 10*" (cgs), indicating ~63% of the immedi- 
ate post-shock momentum has been translated into pressure 
ahead of the contact discontinuity. 

6.5 Evidence for radiation from the shocked gas 

If the shocked outflow scenario does apply to NGC 4051 
the question arises as to whether there is direct evidence for 
radiation from the post-shock cooling? 

The high temperatures in the shock require strong cool- 
ing which is likely to be dominated by Compton scattering 
of the AGN's radiation field. King (2010) finds that this typ- 
ically will have a Compton temperature Tc ~ 10^ K, com- 
pared with the much higher adiabatic shock temperature of 
mpv'^/k ~10" K. 

In assessing whether there is any direct evidence of such 
cooling radiation, we recall that Uttley et al. (2003) identi- 
fied a quasi-constant soft X-ray component, modelled as a 
r~3 power law, in a Chandra TOO observation made dur- 
ing a 6- week low fiux state of NGC 4051. The overall X-ray 
spectrum in the present low flux revl739 data is very simi- 
lar to the Chandra observation, and also to that during the 
XMM-Newton observation in 2002 which followed a 20-day 
low flux period. Furthermore, regular monitoring with Swift 
during the present XMM-Newton campaign indicated that 
revl739 also followed a low flux state lasting for several days. 
We therefore take the revl739 spectrum of NGC 4051 to be 
typical of a possible base level spectrum. 

Figure 11 shows the pn data from revl739 fitted with 
two continuum components. Above the sharp spectral break 
at ~l-2 keV the spectrum is parameterised by a hard power 
law (r~0.8), while the steep soft X-ray component can be 
modelled by either a much softer power law (r~3.5), or a 
Comptonised spectrum (as shown in the figure) with kT~0.3 
keV and optical depth r~0.3. A strong Gaussian emission 
fine (equivalent width 250±30 eV) at 6.38±0.01 keV is con- 
sistent with Fe K fiuorescence from the hard spectral com- 
ponent being reflection dominated, as reported in Po04. 

The short-term flux variability seen throughout the rest 
of the 2009 XMM-Newton observation of NGC 4051 is no- 
tably absent in rev 1739 (Vaughan et al. 2010), supporting 
the view that the soft X-ray component in flgure 11 does in- 
deed comprise a quasi-constant low flux emission, with the 
lack of variability indicating its origin in an extended region. 
The post-shock cooling shell could flt that description and it 
is notable that the present RGS spectra show that absorp- 
tion is limited to self-absorption in the broad emission lines 
at the lowest continuum fluxes. The luminosity in the Comp- 
ton component of flgure 11 ( 0.2 keV) is ~3 x 10''^ erg s~^, 
comparable to the mechanical energy lost in traversing the 
post-shock region. 

We have suggested that the broad emission lines and re- 
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Figure 11. Parametric model fit to the lowflux revl739 pn data. 
The hard power law and soft Comptonised continuum compo- 
nents are shown as dotted and dashed lines 

combination continua (RRC) observed in the RGS spectrum 
arise from additional 2-body cooling of the post-shock out- 
flow. Fitting in XSPEC finds fiuxes of ~1.5 x lO""* photons 
cm~^ s~^ and ~2 x 10~^ photons cm~^ s~^, respectively, 
for the broad emission lines of OVII and OVIII Lyman-a. 
Addition of the corresponding broad emission lines of Ne, N 
and C yields a total soft X-ray fiux of > 10""^ photons cm~^ 
s~^, and a luminosity >2 x 10''° erg s"^. Strong RRC of 
O, N and C and emission from Fe-L lines increase the total 
observed recombination cooling to ~5 x 10*° erg s~^, mak- 
ing a significant addition to the cooling in the later stages 
of the post-shock fiow. 

We conclude that identifying a quasi-constant soft con- 
tinuum and soft X-ray emission features in NGC 4051 with 
the cooling of the shock outflow is reasonable on energetic 
grounds. 

6.6 Relating NGC 4051 to other AGN with 
powerful ionised winds 

Evidence has been growing for ultra-fast winds of highly 
ionised matter in a number of AGN, carrying mechanical en- 
ergy up to 10% of Lto! ( Pounds and Reeves 2009). The direct 
determination of a large-angle flow in the QSO PG1211-I-143 
was important in confirming that the v~0.13c wind was en- 
ergetically significant in terms of AGN feedback. Indirect 
support for such fast outflows to be typically of wide an- 
gle has recently been obtained in a survey of bright AGN 
by Tombesi et al. (2010), who find some 30% of their sam- 
ple show evidence for an ionised wind of v~0.1c. As noted 
above, the potential importance of such powerful winds lies 
in providing a feedback mechanism linking the growth (and 
termination of growth) of supermassive black holes in AGN 
with that of their host galaxy. This remains true if the im- 
pact is actually delivered by a momentum-driven thrust, as 
would be the case if the initial ionised wind lost much of 
its mechanical energy in shocks before reaching the star- 
forming region. The analysis outlined above suggests that 
may be the case for the bright, nearby Seyfert galaxy NGC 
4051. 

The question of how a highly ionised gas is accelerated 
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to such high velocities is clearly important. Continuum driv- 
ing for AGN accreting at a super-Eddington rate was pro- 
posed by King and Pounds (2003), and appears readily ap- 
plicable to PG1211-I-143. However, in the case of many of the 
Tombesi et al. sample, including NGC 4051, current mass 
estimates suggest they are accreting more typically at only 
~10-20 % of LEdd- 

If the Tombesi et al. findings are confirmed, with highly 
ionised outflows at v^O.lc being relatively common in bright 
nearby galaxies, then perhaps Eddington or mildly super- 
Eddington accretion is also more common than generally 
believed. A case for AGN black hole masses being over- 
estimated has rcccntlj' been argued by King (2010a). Fur- 
thermore, if the ejection of fast outflows is intermittent, as 
suggested here for NGC 4051, then only where such a wind 
is current or was launched very recently will it retain a line- 
of-sight column close to the theoretical value for continuum 
driving of Nh~10^* cm'^ 



6.7 Is the forbidden line emission dominated by 
swept up ISM? 

Strong and narrow emission lines in the RGS data, of similar 
strength in both low and high flux spectra, arise from the 
'forbidden' transitions in the ls-2p triplets of He-like NVI, 
OVII and NelX. In each case we find the line is at best only 
marginally resolved, with the strongest OVII line indicat- 
ing a FWHM<250 km s~^. When adjusted for the known 
redshift of NGC 4051 the OVII forbidden line also has a 
very low outfiow velocity of -125±40 km s~^. We note that 
the velocity width is consistent with the [OIII] line width 
in the NLR, of 210-330 km s"^ FWHM (De Robertis and 
Osterbrock 1984). 

In P04 constraints on the low ionisation/low velocity 
gas were obtained from the OVII emission line flux by not- 
ing that the 2002 November XMM-Newton observation took 
place some 20 days after the source entered an extended low 
flux state, while the emission line strength of the OVII for- 
bidden line was essentially the same as when NGC 4051 
was much brighter in 2001 May. This was taken to indi- 
cate a recombination time ^ 2 x 10''s, and a plasma density 
n<8 X 10*cm~^. We now note the forbidden line fluxes in 
the 2009 RGS data are consistent with those measured 7 
years earlier. If that is a real measure of lack of variability it 
would indicate a still lower density n<10'\ iii^'\ or a physi- 
cal extent of the forbidden line emission region larger than 
the present estimate for a post shock region of radius ~0.3 
pc. 

The question then arises, does a substantial fraction of 
the forbidden line emission come from the slower moving 
swept-up ISM ahead of the outer shock? This possibility 
is raised by noting that the resonance absorption line of 
OVII, identifled in the context of a shocked wind with mat- 
ter ahead at the contact discontirmity, has a lowest velocity 
component of 440±60 km s^^, well separated from that of 
the observed OVII forbidden line. Absorption in OV and - 
less unambiguously - in OIV also have lower velocities and 
lie well to the low ionisation side of the linear correlation 
with velocity that fits the main outflow stages in figure 8. 

Estimating the baryon mass of the swept-up ISM firom 
the present data depends both on the assumed density and 



also - quite strongly - on the relevant ionisation parame- 
ter. While the optimum ionisation parameter for OVII is 
log^~0.9, the very low velocity would associate it more with 
absorption from OIV- VI and log^~0, where OVII would rep- 
resent only ~3 % of total oxygen. 

With the measured OVII forbidden line flux of ~1.5 x 
10"* photons cm"^ s"\ a distance to NGC 4051 of 15.2 
Mpc, 3% of oxygen in the form of OVII, and a recombination 
rate at kT -^3 eV of 2 x 10"" cm^ s"^ (Verner and Ferland 
1996) we find an emission measure of '^lO^^cm^'^. 

For a particle density n-^lO^cm"^ that emission mea- 
sure would correspond to a total swept-up mass of ~7500 
Mq. Assuming an ISM density close to the virial equilib- 
rium value (i.e. the isothermal sphere value for cr = 88 km 

in NGC 4051), the swept up mass within r~10^* cm, 
for a gas fraction of 0.15, would be ~3.5 x 10^* gm (1.7 x 10^ 
Mq. The implication is that previous Eddington episodes 
in NGC 4051 have substantially reduced the gas density in 
the inner core of the galaxy. 



7 SUMMARY 

Extended new XMM-Newton observations of the Seyfert 1 
galaxy NGC 4051 have shown a rich X-ray absorption line 
spectrum from 16 K-shell ions from C to Fe. The increased 
sensitivity and energy span of these observations reveal a 
velocity-structured outflow ranging over at least two orders 
of magnitude in velocity and ionisation parameter. We find a 
strong correlation of velocity and ionisation parameter, with 
both parameters being structured in at least 3 regimes. 

We assess the observations in terms of a mass-conserved 
outfiow, where a highly ionised, high velocity wind is slowed 
on impact with the local interstellar medium or slower mov- 
ing ejecta, resulting in a strong shock. While most of the 
initial mechanical energy is lost in the shock, the total out- 
fiow momentum is maintained. 

On this picture, we associate the intermediate velocity 
components, observed at ~3000-7500 km s~^ and with the 
highest opacity at intermediate ionisation parameters, with 
the rapidly cooling shell (or shells) of post-shock gas, while a 
well-separated continuum component at ^400-700 km s~^, 
in a lower ionisation state, is attributed to the gas accumu- 
lating ahead of the contact discontinuity. We speculate that 
a still lower velocity component may be in the ISM swept 
up by the forward shock. 

Compton cooling of the shocked gas by the AGN radi- 
ation field is predicted to produce a soft continuum emis- 
sion containing much of lost outflow kinetic energy and we 
propose that part of this might be observed as a residual, 
quasi-constant soft X-ray component (r~3.5) in the lowest 
flux spectra of NGC 4051. 

The broad line emission and strong RRC observed from 
several of the most abundant ions in the RGS spectrum pro- 
vide additional cooling as the post-shock flow recombines. 
Self-absorption arising in the higher line-of-sight column 
density around the circumference of a thin shell of post- 
shock gas could explain the persistence of low velocity ab- 
sorption, across a wide range of ions, even when the contin- 
uum is faint. 

If a non-varying OVII forbidden line is a measure of the 
swept-up ISM ahead of the forward shock, the estimated 
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mass of ~7500 represents only ~4% of that for a virial 
gas density, suggesting that previous Eddington episodes in 
NGC 4051 have substantially reduced the gas density in the 
central regions of the galaxy. 
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